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ABSTRACT 

We present the first results from the Hubble Space Telescope {HST) Advanced Camera for 
Surveys (ACS) Grism Parallel Survey, a large program obtaining deep, slitless ACS grism spec¬ 
troscopy of high-latitude HST parallel fields. We report on 11 high Galactic latitude fields here, 
each with grism integration times > 12 ks. We identify 601 compact emission line galaxies at 
z < 1.6, reaching emission lines to a flux limit of > 5 x 10“^® ergscm”^ s“^ (Sc)- We determine 
redshifts by cross correlation of the target spectra with template spectra, followed by visual in¬ 
spection. We measure star formation rates from the observed [O II] A3727, [O III] A5007, and Ha 
line fluxes. Follow-up observations with the Keck telescope of one of the survey fields confirms 
our classification and redshifts with a{z) ~ 0.02. This is one of the deepest emission line surveys 
to date, covering a total area of 121 arcmin^. The rough estimate of the co-moving number 
density of emission-line galaxies in our survey at 0.3 < z < 1.3 is ~ 4.5 x 10“® h^Q Mpc“®. We 
reach deeper into the emission-line luminosity function than either the STIS or NICMOS grism 
parallel surveys, finding an apparent space density of emission line galaxies several times higher 
than those surveys. Because of the AGS high spatial resolution, our survey is very sensitive to 
faint, compact galaxies with strong emission lines and weak continua. The AGS grism survey 
provides the co-moving star formation density at z < 1.6 at a high level of completeness. 

Subject headings: galaxies: catalogs — galaxies: high redshift — surveys 


1. Introduction 

One of the key unanswered questions of modern 
cosmology is the origin and extent of the decline 
in the star formation rate of the Universe between 
z ~ 1.5 and the present epoch (e.g., Gallego et al. 
1995; Tresse & Maddox 1998; Lilly et al. 1996; 
Gowie et al. 1999; Hippelein et al. 2003). Based on 

^ Spiteer Science Center, California Institute of Technol¬ 
ogy, MS 100-22, Pasadena, CA 91125 

^Astronomical Institute of St.Petersburg University, 
Petrodvoretz, 198504, Russia 

®Massachusetts Institute of Technology, Cambridge, 
MA 02138 

^Jet Propulsion Laboratory, California Institute of 
Technology, Mail Stop 169-506, Pasadena, CA 91109 

® Steward Observatory, University of Arizona, 933 N. 
Cherry Ave., Tucson, AZ 85721 

^Department of Astronomy, 601 Campbell Hall, Univer¬ 
sity of California, Berkeley, CA 94720 


studies of star-forming and evolving galaxies as a 
function of look-back time, a picture has emerged 
in which more massive galaxies undergo a larger 
fraction of their star formation at earlier times 
than less massive ones (e.g., Gowie et al. 1996; 
Gavazzi & Scodeggio 1996; Cimatti et al. 2004; 
Glazebrook et al. 2004; Kodama et al. 2004). This 
potentially presents a challenge for existing mod¬ 
els of galaxy formation (e.g., Brinchmann et al. 
2004). To fully address this question, it is nec¬ 
essary to obtain a comprehensive sample of star¬ 
forming galaxies of different masses and morpho¬ 
logical types at a broad range of redshifts; deep 
spectroscopic surveys provide ideal samples for 
this endeavor. 

The low sky background and high spatial reso¬ 
lution afforded by space makes emission line sur¬ 
veys from above the atmosphere particularly sen¬ 
sitive and powerful tools for studying galaxy for- 


1 



mation and evolution. Ground-based, objective- 
prism programs, such as the Kitt Peak Inter¬ 
national Spectroscopic Survey (KISS; Gronwall 
et al. 2004) have been quite effective at identi¬ 
fying bright, low-redshift Ha emitters over wide 
areas of sky. To date, KISS has identified 2266 
emission-line objects over 182 deg^, to a limiting 
flux of 1X 10“^^ ergs cm“^ s“^. Ground-based sur¬ 
veys have also been effective at identifying faint, 
high-redshift Lya emitters in gaps between the 
telluric night-sky lines (e.g., Hu et al. 1998; 
Rhoads et al. 2000; Kodaira et al. 2003; Stern 
et al. 2005), reaching typical limiting line fluxes 
of « 2 X 10“^^ ergscm“^s“^. Slitless, grism sur¬ 
veys from space provide the opportunity to iden¬ 
tify fainter lines unobstructed by telluric OH emis¬ 
sion. Furthermore, parallel programs with the 
Hubble Space Telescope (HST) provide, at no cost 
to the observatory efficiency, surveys which are less 
susceptible than pencil beam surveys to the bias 
induced by cosmic variance (e.g., Gohen et al. 
2000 ). 

Previous HST parallel slitless spectroscopic 
programs with the Space Telescope Imaging Spec¬ 
trometer (STIS; Gardner et al. 1998; Teplitz et al. 
2003) and the Near-Infrared Camera and Multi- 
Object Spectrograph (NICMOS; McCarthy et al. 
1999; Yan et al. 1999) have left a valuable scien¬ 
tific legacy studying faint emission line galaxies 
out to high redshift. In particular, the lack of tel¬ 
luric emission lines allowed the NICMOS Parallel 
Survey to identify an impressive census of high- 
redshift galaxies whose optical features are shifted 
into the near-infrared. ACS, though observing 
at the same wavelength regime as STIS, provides 
higher spatial resolution data and is a factor of a 
few more sensitive to emission lines, thus provid¬ 
ing a significant increase in our ability to study 
faint emission-line galaxies (Pirzkal et al. 2004). 

In this paper we present the analysis of eleven 
deep spectroscopic fields from the ACS parallel 
survey. ACS slitless spectra provide unprece¬ 
dented sensitivity in the range 5500 < A < 10500 
A, where ground-based spectroscopy is challenged 
by the night sky. We discuss a sample of z < 1.6 
galaxies selected by the presence of [O II] A3727, 
[O III] A5007 and/or Ha emission features. Our 
survey is well-suited for exploring the faint end of 
the star-forming galaxy luminosity function, iden¬ 
tifying Lya emitting galaxies at 4 ;< z ;< 7, and for 


studying the spatially-resolved star formation rate 
in individual galaxies. This paper presents the 
first results from our survey. In §2 we describe our 
observation and data reduction methodologies. In 
§3 we describe how redshifts are determined and 
present a comparison of our ACS grism redshifts to 
ground-based observations obtained at Keck Ob¬ 
servatory. In §4 we present initial results from our 
survey, including measurements of star formation 
rates for individual galaxies (§4.2) and ACS mor¬ 
phologies of actively star-forming galaxies (§4.3). 
Our results are summarized in §5. The first paper 
in this series describes our data reduction scheme 
in detail (Paper I; Chen et al. 2005); future pa¬ 
pers will address the faint end of the luminos¬ 
ity function for star-forming galaxies (Paper HI; 
Drozdovsky et al., in preparation), and the early- 
type galaxy sample at0.6<z<1.3 detected with 
the grism spectra (Paper IV; Yan et al., in prepa¬ 
ration) . 

Throughout this paper, unless otherwise spec¬ 
ified, magnitudes refer to the Vega system and 
we adopt a flat, A-dominated universe (Hq = 
70/170 km s“^ Mpc“^, Om = 0.3, Oa = 0.7). 

2. HST/ ACS Observations 

All HST data presented here were obtained 
with the Wide Field Camera (WFC) on ACS. 
The observing program was designed to acquire 
a pair of images at each pointing: a direct image 
taken with a broad-band filter (typically F775W 
or F814W) and a dispersed image taken with 
the G800L grism (see Table I). The direct im¬ 
ages are important both for registering the grism 
frames to a common origin and for the 0**' or¬ 
der wavelength calibration of the grism spectra. 
The WFC/G800L grism has a mean dispersion of 
~ 40 A pix ^ in the first order (Pavlovsky 2004). 
The actual, observed resolution is a function of 
the apparent image size convolved with the WFC 
point spread function (PSF). From the observa¬ 
tions described herein, the realized spectral reso¬ 
lution, R\, is ~ 80 to ~ 150. Small variations in 
the PSF due to changes in the optical telescope 
assembly (cf., breathing) and longer term changes 
in the internal structure of ACS introduce small 
variations in the maximal achievable resolution. 
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2.1. Parallel observing mode and schedul¬ 
ing 

Approximately half of the observations for this 
investigation were obtained with HST operating in 
parallel mode as part of GO/PAR program 9468^ 
(P.I. Lin Yan). While this allowed us to collect far 
more data than would be possible in a single pri¬ 
mary program, it limited our ability to plan and 
execute the observations in a manner that opti¬ 
mized the scientific return. Since July 2002, the 
ACS/WFC parallels were scheduled during much 
of the time for which NICMOS or STIS were the 
primary instruments. The observing algorithm for 
this program was quite simple. In each full orbit 
one of the following exposure sequences was se¬ 
lected: F606W and/or F814W imaging followed 
by G800L grism observations with exposure times 
of approximately 500 s. The preferred ratio of ex¬ 
posure times was F814W:F606W:G800L = 1:1:3. 
The selection of the exposure sequence for any 
given orbit was nearly random, but was weighted 
in favor of the spectroscopy. The dither step and 
orientation depended upon the primary observing 
programs, and we preferentially use programs with 
fixed orientation and small dither offsets. 

Another half of our data were obtained from 
the HST Archive as part of guest observer paral¬ 
lel GO/PAR program 9482 (P.I. James Rhoads), 
and data for the J 08:08-1-06:43 field are from GO 
program 9405 (P.I. Andy Fruchter). The basic ob¬ 
serving approaches were similar to ours, though 
with alternate filter selections (see Table 1). 

During the period from July 2002 to Octo¬ 
ber 2003, approximately 800 AGS grism exposures 
(roughly 200 orbits) were observed. These were 
distributed in ~ 40 independent pointings with in¬ 
tegration time of each individual grism exposure 
ranging from 300 s to 1200 s. Rather than adopt a 
strict latitude cutoff we chose to reduce all of the 
data and reject those with high stellar densities. 
Depth of the grism images is the major limitation 
of the study. For this study we choose 11 fields 
with total grism exposure time more than > 12 ks, 
covering about 121 square arcminutes. 

The final depth achieved varied from field to 
field. We define our limiting depths as Scr within 
a 4 pixel aperture. This aperture reasonably repre- 

^ Information about the observations can be gleaned directly 
from the STScI WWW pages linked to the program ID. 
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Fig. 1.— Area-depth histogram for the ten ACS 
grism parallel fields discussed here. Flux limits 
refer to 3cr limits in 4 pix apertures; this aperture 
reasonably represents the ACS/WFC PSF. 

sents the area of the ACS/WFC PSF. Fig. 1 shows 
an area-depth histogram for our survey. Our me¬ 
dian limiting depth is ~ 1.6 x 10“^^ ergscm“^s“^ 
and our deepest three fields reach a depth of 
~ 5 X 10“^® ergs cm“^ s“^, which is approximately 
a factor of two deeper than typical ground-based 
narrow-band/spectroscopic surveys and 2—5 times 
deeper than the HST/NICMOS (McCarthy et al. 
1999; Yan et al. 1999) and STIS (Teplitz et al. 
2003) parallel surveys. 

2.2. Data processing and analysis 

The extraction of the spectra from the grism 
images is a multi-stage process. The main steps 
of the reduction procedure are: standard CALACS 
pipeline reduction (e.g., bias subtraction); cosmic 
ray detection and global sky background removal; 
combining the 2D images; generating object cat¬ 
alogs from the direct images; extraction and cali¬ 
bration of ID spectra from the coadded grism im¬ 
ages; combination of the ID spectra for objects 
observed at different positions and/or with differ¬ 
ent orientations. The software used for the spec¬ 
tral extraction is written by one of us (Chen et al. 
2005, Paper I). Paper I presents the design and 
performance of our data reduction software. We 
detail our reduction steps below. 
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2.2.1. CALACS pipeline and image combination 

We begin with the output from the STScI data 
pipeline, CALACS, which does basic processing of 
the individual 2D frames. All data have first or¬ 
der bias subtraction, dark subtraction, and bad 
pixel masking applied. Direct images are then 
corrected for the flat-held response. Since the flat- 
field response is a function of wavelength for each 
pixel and so depends on the location of sources in 
the dispersed data, the grism data are not flat- 
fielded.^ 

We combine the direct images using a modi¬ 
fied version of the Multidrizzle package (Koeke- 
moer et al. 2002). Multidrizzle is an implemen¬ 
tation of the Blot/Drizzle technique (Fruchter & 
Hook 2002) which provides an automated method 
for distortion-correcting and combining dithered 
images. Multidrizzle corrects for gain and bias 
offsets between WFC chips and identifies and re¬ 
moves cosmic rays and cosmetic defects. The qual¬ 
ity of the image combination relies on accurately 
determining the offsets between images. The orig¬ 
inal program aligns images using the World Co¬ 
ordinate System (WCS) header keywords. Unfor¬ 
tunately, there are times when those WCS val¬ 
ues are inaccurate, resulting in a misalignment of 
the final output drizzled product. We have there¬ 
fore enhanced Multidrizzle with the ability to 
internally verify and update WCS-calculated off¬ 
sets using image cross-correlation or point-source 
matching algorithms. The derived shifts and ro¬ 
tations are determined for the geometrically cor¬ 
rected frames resampled onto a common WCS 
frame. The absolute astrometry of the final driz¬ 
zled images were verified using guide stars as well 
as extragalactic sources from NED. The astromet¬ 
ric uncertainties are dominated by the accuracy 
of the coordinates of the guide stars located in a 
field area and differ from field to field with a me¬ 
dian value of ~ 0"2. 


^The aXe ACS grism reduction code of Pirzkal et al. (2002) 
uses a data cube constructed from observations through 
narrow-band filters, interpolated to the wavelengths of ex¬ 
tracted pixels, to fiat-field grism spectra. This technique 
reduces pixel-to-pixel scatter within the 6000 — 8500 A 
wavelength regime from 3% to 1%, but shows no improve¬ 
ment beyond 8500 A due to a lack of long-wavelength ACS 
narrow-band filters. Since the G800L grism mode remains 
sensitive to Si 10 000 A, we opt to omit fiat-fielding during 
spectral extraction (see Paper I for details). 


For the grism data, our ability to stack raw ACS 
frames is limited by the substantial geometric dis¬ 
tortion and the extent of the offsets between in¬ 
dividual exposures. According to the ACS Instru¬ 
ment Handbook, the plate scale changes by 8% be¬ 
tween the two diagonal corners of the field. This 
discrepancy is less than a pixel for dither offsets 
smaller than 10 pixels. Therefore, we register dis¬ 
torted ACS grism images that have pointing off¬ 
sets less than 10 pixels into substacks, taking care 
to mask bad pixels and applying a 5a clipping al¬ 
gorithm to remove hot pixels and residual cosmic 
rays. For fields with larger offsets or varying orien¬ 
tations between exposures, we create multiple sub¬ 
stacks. Spectral extraction is usually performed 
on these substacks, then the final coadd is done 
for 1-dimensional spectra. 

Our spectral extraction software (Paper I) re¬ 
quires a pair of aligned direct and dispersed im¬ 
ages in the original post-pipeline format; i.e., each 
WFC chip image must reside in a separate FITS 
files with the original ACS distortion. We make 
use of the Blot program, which transforms direct 
images back to the reference frame of G800L im¬ 
age stack(s). We extract ID spectra from the sub¬ 
stacks. In cases with multiple substacks, the ID 
extracted spectra are resampled onto a common 
grid and weight-averaged to create the final deep 
spectrum. 

2.2.2. The Catalog 

We generate object catalogs from the stacked 
direct image using SExtractor (Bertin & Arnouts 
1996), applying a 1.5cr per pixel detection thresh¬ 
old and 8 pixel minimum area requirement. Ta¬ 
ble 2 presents an example of the catalog produced 
from our imaging dataset; the full catalog of emis¬ 
sion line galaxies is produced in the on-line version 
of this paper. Below we detail some of the listed 
parameters: 

• Equatorial coordinates are the barycenter 
position of a source, that were measured in 
a WCS-corrected direct image (see § 2.2.1). 

• Magnitudes are measured in Kron-like AUTO 
magnitudes and transformed to the Vega 
photometric system using following zero- 
points: 26.398 for F606W, 25.256 for F775W, 
25.501 for F814W, and 24.326 for F850LP 
(Pavlovsky 2004). 
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• a,b are the major and minor axis lengths 
of the fitted ellipse (in units of the 0'/05 
ACS/WFC pixels). 

• 6 is the position angle of the major axis with 
respect to the dispersion direction. 

• Concentration, the Abraham et al. (1996) 
concentration parameter, is the ratio be¬ 
tween the flux in the central 30% of the pix¬ 
els as compared to the total flux. 

• Asymmetry, the Abraham et al. (1996) 
point-asymmetry index, is the absolute sum 
of the differences between point symmetric 
pixels divided by the total flux. 

2.2.3. 1-D spectral extraction and calibration 

The extraction and calibration of the ID spec¬ 
tra was carried out with custom software devel¬ 
oped by Hsiao-Wen Chen and described in detail 
in Paper I. We briefly describe the methodology 
here; the interested reader is directed to Paper I 
for details. 

We first align sources in the direct image (po¬ 
sition Xi,yi) with the corresponding position of 
their zeroth orders in the dispersed image (x?, y^), 
adopting the empirically-derived transformation: 

= -122.1591 -b 1.00442 • Xi - 0.00395 • yf, 
y° = 3.6392 - 0.00014 • Xi + 1.00088 • y^ 

( 1 ) 

for WFC Chip 1, and 

a;° = -115.2942 -b 1.00444 • x, - 0.00352 • yp, 
y° = 2.5025 - 0.00028 • Xi + 1.00068 • y^ 

( 2 ) 

for WFC Chip 2. The extraction region of the 
spectrum was then computed using transforma¬ 
tions derived as part of the calibration program. 
Using the object brightness distribution inside its 
segmentation map on the direct image, the cor¬ 
responding 2D spectra were modeled and over¬ 
lapping spectra were de-blended via iterative, 
multiple-profile fitting procedure for all spectra 
in a frame simultaneously. For objects close to 
the edges of the dispersed images, only partial 
spectra were extracted. For faint sources, only 
first order light is considered; for bright sources, 
we improve the signal by including higher orders. 
The rms accuracy of the wavelength calibration for 


G800L is approximately 7A (Pasquali et al. 2002). 
We flux calibrate the extracted spectra using cal¬ 
ibration curves derived by Pirzkal et al. (2002) 
from observations of white dwarfs and Wolf-Rayet 
stars. The accuracy of the spectrophotometry is 
limited by uncertainties in the wavelength cali¬ 
bration, the various detector flat-held effects, the 
object deblending, and variations in the quantum 
efficiency within individual pixels. We estimate 
that the absolute flux calibration is accurate to 
approximately 5% from 5000 A to 9000 A. 

Spectra were also extracted in a parallel ef¬ 
fort using the aXe software developed at ST-ECF 
(Pirzkal et al. 2002). The aXe software follows a 
similar strategy to our code, except it does not 
perform any deblending due to higher orders from 
bright objects. Spectra are extracted in weighted 
boxes, with flat-fielding performed based on ob¬ 
servations through narrow-band Alters, interpo¬ 
lated to the wavelengths of the extracted pixels. 
Flux calibration uses the same calibration curves 
as used by us. The overall results of the extrac¬ 
tion and analysis using the aXe and our software 
were similar, except in cases of overlapping spec¬ 
tra where the de-blending included in our code 
produces significantly cleaner extractions. 

3. Redshift Identifications 

We next cross-correlate the extracted, ID spec¬ 
tra with a set of stellar and galactic spectral tem¬ 
plates in order to deduce an approximate spec¬ 
tral type and redshift for each source. Sources 
were divided into four groups: Galactic stars, 
emission-line galaxies (ELGs), early-type galaxies, 
and galaxies with a Lya break at 4 < 2 < 7. In 
this paper we present analysis of only the ELGs. 
Emission-line fluxes and equivalent widths were 
measured using Gaussian fits to the line and poly¬ 
nomial fitting of the underlying continuum, per¬ 
formed using the ESO-MIDAS Alice package. 
Fig. 2 presents direct images and spectra of a sam¬ 
ple of 6 ELGs from this survey. Galaxies were se¬ 
lected to illustrate the range of data quality. Ta¬ 
bles 2 and 3 contain the catalog of imaging and 
spectroscopic properties of the ELGs detected in 
our fields^. 


®The complete version of these tables are in the electronic 
edition of the Journal. The printed edition contains only a 
sample. 
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Fig. 2.— Typical emission-line galaxies detected from the ACS Grism Parallel Survey, showing direct images, 
dispersed images, and extracted ID spectra. All spectra are shown over observed wavelength range in flux 
units of erg s“^ cm“^ . The dotted lines indicate the best-fitting spectral template. Data are sorted by the 
reliability of their line identifications: the top row shows example of robust, class ‘a’ redshift identifications, 
the middle row shows examples of likely, class ‘b’ redshift identifications, and the bottom row shows less 
secure, class ‘c’ redshift identifications. The postage stamp images are 5" on a side, oriented as per the 
original data. 
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3.1. The identification of emission lines 

Although ~ 80% of objects with visible emis¬ 
sion lines were successfully classified with the cross 
correlation technique, the accuracy of the red- 
shift determination was low and required visual 
verification. We compared the ACS spectra to 
star-forming galaxy templates from Kinney et al. 
(1996) smoothed to the grism resolution. To con¬ 
vince ourselves that any detected emission line 
was real, we eliminated all possible false signals 
(e.g., zero order images, persistent images, cosmic 
rays, etc.). For final spectral classification, three 
criteria were used: identified emission line(s), mor¬ 
phology of the emission line(s), and the continuum 
spectral energy distribution (SED). 

To be specific, three major emission features 
could be present at our observed wavelength 
range: Ha blended with [N II] AA6548,6584 (for 
galaxies at 2 ; < 0.5), the unresolved [O III] 
AA4959,5007 doublet -h H/3 (0.5 < z < 1.0), and 
[O II] A3727 (0.6 < z < 1.7)^. Many sources 
show multiple emission lines, allowing unambigu¬ 
ous redshift identifications. We note that since 
the ratio of [O III] AA4959,5007 flux to both Ha 
and [O II] A3727 flux varies significantly in both 
the local universe (Terlevich et al. 1991; Kenni- 
cutt 1992; Izotov et al. 1994) and at high redshift 
(e.g., Kobulnicky et al. 1999; Teplitz et al. 2000; 
Pettini et al. 2001), line flux ratios do not provide 
a strong redshift indicator. 

Galaxies with only a single emission line require 
more attention. We assume that isolated emis¬ 
sion lines are either Ha or [O II] A3727. Since the 
[O n]:[0 III] ratio can vary from 0.1 to 10, and the 
[O III]:Ha ratio can vary from 0.33 to 1 (Kennicutt 
1992), the non-detection of [O III] in an Ha- or 
[O II]-emitting galaxy is possible. We choose be¬ 
tween Ha and [O II] based on the relative position 
and strength of the line, the continuum SED, as 
well as the brightness and morphology of the host 
galaxy in the direct image. For example, a single 
emission line at wavelengths shortward of 7800 A 
that would imply either an Ha-emitter at a low 
redshift (z < 0.2) or an [O II] emitter at moderate 
redshift (z < 1). An Ha-emitter might be ex¬ 
pected to appear noticeably larger and brighter in 
the direct image, possibly showing resolved struc- 

^ Galaxies with Lyo emission at 2 > 4 are not considered 
here. 


ture. Since the space density of low-redshift, star¬ 
forming compact dwarf galaxies is low, a single, 
red emission line in a faint galaxy (with a blue 
rising continuum at shorter wavelength) is most 
likely [O II] unless there is evidence to the con¬ 
trary. 

Line morphology can also provide a useful red¬ 
shift diagnostic. At our spectral resolution, the 
[O III] line is generally asymmetric due to blend¬ 
ing with H/3, allowing us to secure line identifica¬ 
tions based on the line morphology. On the other 
hand, [N II] AA6548,6584 are not resolved from Ha 
at our resolution, particularly at the low redshifts 
at which these lines are observed. Some of the 
compact low-redshift sources, however, reveal the 
noticeable [S II] AA67I6,673I doublet. 

We assign a quality (reliability) flag to all red¬ 
shift estimations based on the number of detected 
emission lines and significance of their identifica¬ 
tions. Quality flag ‘a’ indicates that there are 
two or more emission lines identified and that 
the [O III]-|-H/3 blend shows a clear asymmetry; 
this flag indicates that the redshift is secure. A 
value of ‘b’ is assigned to galaxies for which there 
is a strong reason for the assignment. Specifi¬ 
cally, such a quality assignment implies we ob¬ 
serve multiple emission lines but doubts remain as 
to their identifications. Quality flag ‘c’ indicates 
even greater uncertainty in the redshift identifi¬ 
cation, typically indicating that only one line has 
been detected. A non-detection of the second line 
can be due to several reasons: the most common 
circumstance is that the galaxy is either at z < 0.2 
or at z > 0.8, and [O III] is in a region of poor spec¬ 
troscopic sensitivity (A < 6000 A or A > 9000 A). 
In a few cases, the spectral range is truncated since 
the galaxy lies near the edge of the field of view. 

3.2. Comparison with Keck spectroscopy 

To test the accuracy of the redshifts and 
study potentially-interesting faint sources, on 
UT 2004 March 19 we obtained spectroscopy 
of two slitmasks targeting ACS-selected galax¬ 
ies (field J 10:03-1-29:06) with the Low Resolu¬ 
tion Imager and Spectrograph (LRIS; Oke et al. 
1995) on the Keck I telescope. These observa¬ 
tions were obtained in non-photometric condi¬ 
tions and integration times totaled one hour per 
slitmask, split into three dithered 1200 s expo¬ 
sures. LRIS is a dual-beam spectrograph: we 
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Fig. 3.— Keck/LRIS spectrum (thin line) of one 
of the ELGs discovered in this survey, overlaid 
on the ACS grism discovery spectrum (thick grey 
line). Both LRIS and ACS spectra produced sim¬ 
ilar redshift determinations, within the expected 
range considering the ~ 25 times lower spectral 
resolution of the ACS grism data. 

used the D680 beam-splitter, the 300 lines mm“^ 
grism (Ablaze = 5000 A; AAfwhm = 9.0 A) on 
the blue arm, and the 400 lines mm“^ grating 
(Ablaze = 8500 A; AAfwhm = 6.4 A) on the red 
arm. Data were processed using standard tech¬ 
niques. Because the night was not photometric, 
we used archival sensitivity functions dating from 
March 2002 for relative flux calibration. 

Unfortunately due to poor weather conditions, 
we were unable to study the more extreme, faint 
sources identified from ACS. However, the LRIS 
data proved quite useful for verifying the accuracy 
of the ACS-derived redshifts. A total of 11 ACS- 
selected ELCs were observed, out of which one 
was of quality ‘a’, five were of quality ‘b’, and five 
were of quality ‘c’. In all cases, the ACS and Keck 
redshifts were consistent. The mean redshift dif¬ 
ference is negligible, (zACS— -^Keck) = — 0.01±0.02. 
Fig. 3 presents both the ACS and the Keck spectra 
of one of the sources, J 10:03:51.318-1-29:05:00.41, 
at redshift z = 0.337. Note that the H/3/[0 III] 
complex is barely resolved by ACS, but is well- 
resolved by Keck. 

Artificial object tests and comparison of our 
redshift estimations with ones measured using the 
higher spectral resolution Keck/LRIS spectra in¬ 
dicate that manual inspection correctly finds the 
redshifts for over 95% of the spectra. 



Log [N II]6584/Ha 

Fig. 4.— The location of three ELCs (diamonds) 
based on the Keck/LRIS spectra in the BPT line- 
ratio diagram from Brinchmann et al. (2004). Two 
lines shows the empiric classification of the ELCs 
by activity type. 

The Keck/LRIS spectra confirms the high effi¬ 
ciency of the ACS grism selection method, which 
yields a high fraction of strong-lined galaxies. 
Namely, 7 out of 11 observed ELGs (64%) have 
EW([0 III] A5007)> 100 A. Such galaxies tend 
to be either star-burst galaxies with little extinc¬ 
tion, or AGN with high-excitation spectra. The 
general way to classify the ELGs by activity type 
(AGN versus starburst) is based on the emission 
line ratios (Baldwin et al. 1981, hereafter BPT). 
We classify three of these strong-lined sources 
as star-forming galaxies based on their [O III] 
A5007/H/3 versus N IIA6584/Ha line-ratio, shown 
in Fig. 4. For the remaining galaxies, compar¬ 
ison of their [O III] A5007/H/3 ratio with their 
[O II] A3727/[0 III] A5007 ratio, immeasurably 
weak [Ne III] A3826 line, and narrow emission lines 
confirm that they are also starburst systems. 

4. Results 

4.1. Areal coverage and depth 

Based on the analysis of eleven fields, we de¬ 
tect 601 galaxies revealing significant emission 
features, corresponding to a surface density of 
~ 5 ELGs arcmin“^. Among them we identify 
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Fig. 5.— Redshift distribution (top) and ob¬ 
served equivalent width distribution (bottom) for 
detected ELGs. 

166 galaxies with Ha emission, 406 galaxies with 
[O III]-|-H/3 emission, and 401 galaxies with [O II] 
emission. Our survey is most sensitive to emission 
line sources with angular sizes (in the dispersion 
direction) of less than 1 arcsec and with broad¬ 
band magnitudes of F814W < 26 mag. Fig. 5 
shows the redshift and equivalent width distribu¬ 
tion for the survey and Fig. 6 plots line flux fluxes 
against F814W magnitude. Our ELG sample has 
a median redshift of 0.66. At low redshifts, the 
fall off in Ha-emitting galaxies is attributed to 
the small volume covered by survey. At high red- 
shifts, the number of sources falls as the sensitivity 
drops and our primary rest-frame optical features 
enter the near-infrared. There is also some plunge 
in the number of detected single-line galaxies at 
0.4 < z < 0.5 and 1.0 < z < 1.1, when second Ha 
and [O III] lines move into the near-IR. 



36 34 33 30 18 

F814W (Vega mag) 


Fig. 6.— Emission line fluxes versus F814W mag¬ 
nitude. The symbols indicate different emission 
features: asterisks indicate Ha emission, empty 
squares indicate ]0 III] emission, and filled circles 
indicate [O II] emission. The ACS grism survey 
identifies galaxies to very faint continuum bright¬ 
ness levels. 

4.2. The emission-line luminosity 

The distribution of the emission-line luminosi¬ 
ties, ignoring dust extinction corrections, is pre¬ 
sented in Fig. 7. The median Ha line luminosity 
of galaxies in our survey is 2.7 x 10"^° h^Q ergs s“^; 
this is 26 times fainter than characteristic luminos¬ 
ity = 7.1 X lO"*^ h^Q ergs s“^ derived from sur¬ 
veys of the local universe (e.g., Gallego et al. 1995; 
Tresse & Maddox 1998), and 132 times fainter 
than = 3.6 x 10^^ h^Q ergs s“^ at z = 1.3 
(Yan et al. 1999). The median (O II] line luminos¬ 
ity is 8.5 X I0'^°/ifg^ ergs s“^. The median ]0 II] 
luminosity for the z > 0.6 sample is two times 
fainter than the local ^[Oii] ~ ^ 10'^^ h^Q ergs 

s“^ derived by the Universidad Complutense de 
Madrid survey (UCM; Gallego et al. 2002). The 
median [O III]-|-H/3 line-blend luminosity can be 
compared with single ]0 III] A5007 line lumi¬ 
nosity Lpjjjj 5 pp.j. = 1 X lO"*^ h^Q ergs s“^ calcu¬ 
lated by Hippelein et al. (2003) for galaxies in 
the redshift range 0.62 < z < 0.65. Our cu¬ 
mulative [O III]-l-H/3 value of = 7.4 x 

10^^° h^Q ergs s“^ is again 14 times fainter. Our 
survey is clearly reaching down to the faint end of 
the emission-line luminosity functions. Our sam¬ 
ple selects star-forming galaxies over a wide range 
of luminosity, from faint emission-line galaxies at 
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Fig. 7.— Histogram of Ha (solid line), [O III]+H/3 
(dotted line) and [O II] (dashed line) luminosities 
in the ACS Grism Parallel Survey. Vertical lines 
indicate the characteristic Ha line luminosity for 
the local universe (Gallego et al.l995) and at z = 
1.3 (Yan et al.l999). 

low redshifts to luminous L* galaxies at high red- 
shifts. 

4.3. The emission-line luminosity and star 
formation rate 

Ha emission is a classic indicator of star for¬ 
mation because it traces the ionizing flux from 
hot stars. Assuming Case B recombination and a 
Salpeter initial mass function over the mass range 
0.1 < M/Mq < 100, we adopt the calibrations 
of Kennicutt (1998) between star formation rate 
(SFR) and Ha luminosity: 

SFRhq, (At© yr“^) = 7.9 x 10“"^^ • Lna (ergs“^). 

( 3 ) 

Assuming the average Ha to [O II] flux ratio of 
0.45, the Kennicutt (1998) relation implies 

SFR[oii] (Afeyr"^) = 1.4x IO^'^A^jqjjj (ergs"^), 

( 4 ) 

subject to the considerable scatter in this flux ratio 
found in surveys of local galaxies, likely associated 
with variations in the metallicity and star forma¬ 
tion histories of individual galaxies. Furthermore, 
we note that since no extinction corrections have 
been applied, the derived luminosities and star for¬ 
mation rates should be considered lower limits. 

While H/3 line is also a good tracer of star for¬ 
mation in galaxies (e.g., Kennicutt 1983), it is 
blended with the (O III] AA4959,5007 doublet in 



log L(Ha) (erg/s) log L[0I1] (erg/s) 


Fig. 8.— [O III]-|-H/3 line-blend luminosities plot¬ 
ted against Ha and [O II] luminosities for the ACS 
Grism Parallel Survey. Lines show the best fit cor¬ 
relation derived from orthogonal regression. While 
the scatter is considerable, a correlation is obvious. 

our low-resolution spectra. The strength of the 
oxygen lines, however, also correlate with the for¬ 
mation rate of massive stars, and may be used as 
a proxy for the SFR in cases where [O II] and Ha 
are unavailable (e.g., Teplitz et al. 2000; Hippelein 
et al. 2003). Due to its high ionization level, the 
luminosity of [O III] depends strongly upon the 
temperature of the ionized gas, which in turn de¬ 
pends upon the metallicity of the galaxy. In order 
to convert the [O III]-|-H/3 line blend into a star for¬ 
mation rate, an averaged intensity ratio between 
these lines and Ha must therefore be established. 
Having detected and measured a large number of 
ELGs, revealing either the [O III] and Ha or the 
]0 III] and (O II] line pairs, we can test the cor¬ 
relation between [O III] luminosity and SFR with 
the understanding that the analysis is inherently 
limited by the unknown properties of the galaxies. 

Fig. 8 compares the [O II], [O III], and Ha lu¬ 
minosities. Assuming the average Ha:(0 II] ratio 
of ~ 0.6, the median ]0 III]:Ha ratio of w 1.3 
we derive from 127 ACS grism ELGs at z < 0.6 
is in agreement with the median ]0 III]:[0 II] 

« 2.2 ratio we derive from 245 ACS grism ELGs 
at0.5<z<1.0. Using these ratios and equations 
(3) and (4), we derive 

SFR[oiii]+h/3 (Al© yr ^)~6xl0 "‘^•L[oiii]-fH/3 (ergs ^). 

( 5 ) 

We expect this SFR indicator to have the largest 
scatter of the three considered; consequently, we 
only use it for the 47 galaxies in our sample (less 
than 8%) which reveal only a single [O III] line 
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Fig. 9.— SFR derived from Ha (circles), [O III] 
(squares) and [O II] luminosities (asterisks). Size 
of the symbol is proportional to the object 
absolute broad-band magnitude (in F814W or 
F775W). The current Ha-derived SFR for the 
Large Magellanic Cloud (LMC) is indicated by an 
open diamond (Kennicutt, Tamblyn & Congdon 
1994). A filled square illustrates the [O II]-derived 
SFR of the dusty, red starburst galaxy HRIO (Dey 
et al. 1999). 

and lack other emission lines. 

We plot the derived star formation rates® 
against redshift in Fig. 9. At each redshift, our dy¬ 
namic range in observed SFR is approximately 1.5 
orders of magnitude, and we see the expected bias 
of finding higher SFRs at higher redshifts. The 
majority of the surveyed Ha emitters are mildly 
star-forming, local galaxies. The median SFR of 
Ha-emitters in our survey is 0.2A^Qyr“^. After 
correcting for [N II] AA6548,6584 contamination 
using the average [N II]:Ha ratio of 0.3 derived by 
Gallego et al. (1997) from a local sample of galax¬ 
ies, this rate is reduced to ~ O.lA 40 yr“^. This 
SFR is typical of local dwarf galaxies: e.g., the 
median SFRs of the local blue compact galaxies is 
about O. 3 A 40 yr“^ (Hopkins et al. 2002), and the 
Ha-derived SFR for the nearest irregular galax¬ 
ies, LMC and SMC are, respectively, 0.26 and 
0.046 A 40 yr“^ (Kennicutt et al. 1994). 

Using equation (4), we determine that the SFR 
of [O II]-emitters in our survey span nearly 3 or¬ 
ders of magnitude, from approximately a few times 

®Note that since our fluxes are not corrected for any absorp- 
tion effect, the values we present are a low limit to the true 
SFRs. 


10 “^ to several times 10 ^ Mq yr”^, with a median 
SFR of about 1 . 2 Af 0 yr“^. Our survey is sensi¬ 
tive enough to detect objects with SFRs as low as 
lA 40 yr“^ up to z « 1 . 2 . At higher-redshift the 
ACS grism survey reaches typical starburst galax¬ 
ies with SFRs of 20 — 6 OA 40 yr“^ (e.g-. Glaze- 
brook et al. 1999; Savaglio et al. 2004). Even 
a distant analog of the ULIRG population, such 
as ERO J164502-I-4626.4 — also known as HRIO 
(SFR[oii] ~ lOAl 0 yr“^; Dey et al. 1999) — can 
be among our high-redshift sample. 

A major concern remains the large uncertainty 
in SFRs derived from oxygen fluxes. It has pre¬ 
viously been noted in surveys of the local uni¬ 
verse that the [O II]:Ha ratio correlates with to¬ 
tal galaxy luminosity, such that brighter galaxies 
have lower [O II]:Ha ratios (e.g. Jansen et al. 2001; 
Tresse et al. 2002). This correlation is thought 
to be related to metallicity. The variation in the 
[O III]:Ha ratio is even more dependent on metal¬ 
licity, as well as on the effective temperature of 
the gas and the ionization parameter (Kennicutt 
et al. 2000 ), so we proceed with this caution in 
mind. We note, however, that the [O II] and 
[O III] measurements of the SFR do not show a 
large discontinuity with the Ha measurements at 
the transitional region of 2 ~ 0.5 (Fig. 9); appar¬ 
ently, the SFR determinations from oxygen lines 
are not completely wrong. 

A rough estimate of the number density of 
star-forming galaxies can be made using galax¬ 
ies detected at 0.3 < z < 1.3, such that 
their Ha and [O II] lines are in efficient re¬ 
gions of the sensitivity curve. The total an¬ 
gular area of the eleven selected fields is 121 
arcmin^, corresponding to a co-moving volume 
of 1.13 X 10® hjQ Mpc®. We detected 506 ELGs in 
this redshift interval, giving a co-moving number 
density of 4.5 x 10“® Mpc“®. The co-moving 
number density of [O II]-emitters in our survey 
at 0.5 < z < 1.3 is ~ 3.8 x 10“®/i^g® Mpc“®. 
This density is about 7 times higher than that de¬ 
tected by the STIS parallels (Teplitz et al. 2003) 
at 0.5 < z < 1.2. The NIGMOS parallel survey 
(McGarthy et al. 1999) detected 33 Ha emitters 
at 0.7 < z < 1.9 in the co-moving volume of 
0.78 X 10®h7g® Mpc®. Thus the co-moving num¬ 
ber density of ELGs averaged over their volume 
is ~ 0.4 X 10“®/lyg® Mpc“®. As Yan et al. (1999) 
and Teplitz et al. (2003) point out, their surveys 
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Fig. 10.— Half-light radius versus total magni¬ 
tude for three ELG redshift bins (large symbols, 
as indicated). For comparison, all sources in our 
direct images are plotted as small dots. The con¬ 
centration of the ELGs among the objects of small 
angular size objects is evident. 

probed starburst galaxies from the upper end of 
the luminosity function. Our data set is able 
to detect much fainter emission lines and measure 
the co-moving star formation density at z < 1.6 to 
better completeness levels. On the other hand, the 
redshifts measurements based on the ACS/WFG 
slitless spectra are subject to larger uncertainties 
due to the lower spectral resolution of ACS as 
compared to STIS and NICMOS®. 

4.4. The morphology of the ELGs 

The high spatial resolution of the ACS/WFG 
images enable a detailed study of the morpho¬ 
logical properties of ELGs. We defer a detailed 
study to a future paper and present here plots 
of basic morphological parameters for the ELG 
sample, such as half-light radius, concentration, 

®The deep, targeted NICMOS survey of the 4.4 arcmin^ 
Groth-Westphal strip (Hopkins et al. 2000) detected 37 Ha 
emitters at 0.7 < z < 1.8, corresponding to the co-moving 
number density of ~ 7.3x 10“® Mpc“®. We also expect 
that the GRAPES AGS grism survey of the Hubble Ultra 
Deep Field (Pirzkal et al. 2004), with total grism time ~ 28 
hours, will reveal a much larger number density of emission¬ 
line galaxies within the same redshift range as our survey. 



Fig. 11.— The distribution of ELGs on the cen¬ 
tral concentration — asymmetry (C — A) plane, as 
measured in F775W/F814W images. Large sym¬ 
bols represent ELGs; dots show all sources de¬ 
tected in the direct images. The dotted line is 
a rough border between disk and bulge dominated 
galaxies, defined from the visual classihcation of 
well-resolved, bright galaxies. 

and asymmetry. We analyze F775W/F814W im¬ 
ages (observed /-band), which correspond to rest- 
frame B-band for the median redshift of our ELG 
sample, (z) ~ 0.7. The morphology at this rest- 
frame wavelength is dominated by star-formation 
regions. Eig. 10 plots half-light radius against to¬ 
tal magnitude for all objects detected in our imag¬ 
ing, with identihed ELGs indicated. Because of 
band-shifting effects and surface brightness dim¬ 
ming, plot symbols are keyed to redshift. As ex¬ 
pected, the survey is biased towards obtaining red- 
shifts for compact galaxies. 

The morphological classifications of the ELGs 
galaxies were performed on the basis of both the 
visual inspection and automatic classification us¬ 
ing central concentration (C) and rotational asym¬ 
metry (A) indexes (Abraham et al. 1994; Abra¬ 
ham et al. 1996). As shown by Abraham et ah, 
these parameters are remarkably robust to image 
degradation resulting from increased line-of-sight 
distance. Visual classification, however, is more 
sensitive to peculiar and merging galaxies. The 
comparison of visual with quantitative morpholog- 
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ical classifications show good agreement in ~ 70% 
of galaxies brighter than F814W= 22. For fainter 
galaxies, it is not trivial to distinguish E from SO 
galaxies, and Sb or Sc spirals from Sd and Irr 
galaxies. Therefore, based on the visual calibra¬ 
tion of well-resolved, bright galaxies, we distin¬ 
guish here two classes of objects: disk-dominated 
and bulge-dominated galaxies. 

The distribution of the entire imaging data set 
on the C — A plane is shown in Fig. 11, reveal¬ 
ing the variety of their morphological types. The 
fraction of disk versus bulge dominated objects 
in our sample is almost equal. The C — A plot 
is also characterized by a correlated distribution 
of the ELGs, which can be well approximated by 
log% = (0.53 ± 0.03) • logC - (0.08 ± 0.02). The 
small fraction of ELGs among highly asymmetric 
galaxies might be explained by our bias towards 
galaxies of small angular size. The less concen¬ 
trated disk-like systems show a larger spread in 
asymmetry, with many of them resembling inter¬ 
acting, merging, and peculiar systems. 

4.5. ACS simulations and completeness 
tests 

There are several aspects of the data which neg¬ 
atively impact our ability to identify genuine spec¬ 
tral features. Since we are observing in slitless 
mode, each object produces zero, first and second 
order spectra. The first order spectra contain the 
most useful data; the second order is profitable for 
bright objects; zero, negative and higher (> 2"'^) 
orders are sonrces of confusion. The zero-order 
images can can be mistaken for emission features, 
particularly when they fall on the first order con- 
tinua of other objects. The zero-order images are 
slightly dispersed (~ 650 A pix“^) and often ap¬ 
pear bimodal for point sources. The displacement 
between the zero order images and the center of 
the corresponding direct image (near the start of 
the first order spectrum) is ~ 5'.'5, or ~ 110 pix, 
and generally can be identified by matching them 
with either the first order spectra or with an ob¬ 
ject in the direct image. There is a small portion 
of the detector, however, for which zero-order im¬ 
ages can appear without either first order spectra 
or images in the direct frame. In addition to the 
confusion caused by spectral orders, there are ar¬ 
tifacts associated with defective pixels and cosmic 
rays. 


We performed completeness tests using artifi¬ 
cial object trials with the SLIM program (Pirzkal 
et al. 2001), using a set of template starburst 
(SBl) galaxy spectra from Kinney et al. (1996). 
The original SLIM configuration parameters were 
adjusted to correspond to the current best-known 
AGS/WFC/G800L grism spectral trace and dis¬ 
persion descriptions. A total of ~ 600 pairs of 
the synthetic direct images and 2 = 0.7 SBl 
template spectra were generated with SLIM and 
randomly added to our F814W and G800L data 
for the J 13:58-1-62:39 field. The synthetic im¬ 
ages used Gaussian brightness profiles with radii 
071 < r < 075 and magnitudes in the range 23.0 < 
F814W < 24.5. Object spectra were extracted us¬ 
ing the same method we used for the real objects 
and the emission-line parameters were measured 
and compared with input ones. 

Objects were considered recovered if they were 
found in both direct and grism images, the emis¬ 
sion lines were detected at the > 3 (t level, and 
the measured equivalent widths did not exceed the 
initial, injected values by 20 A. Fig. 12 shows the 
outcome of these completeness tests. These plots 
reveal strong selection effects are present in our 
sample of ELGs. More detailed tests, necessary 
for the calculation of the luminosity function and 
measurements of co-moving star formation rate, 
will be discussed in Paper III. 

The quality of the flux and wavelength cali¬ 
bration have also been tested on data taken at 
different epochs and located at different parts of 
the AGS/WFC detector: we find good agreement. 
The rms accuracy of redshifts derived from repeat 
observations is ~ 0.02. 

4.6. Notes on Individual Fields 

16.1. Field J 13:58+62:39 

This field has the largest total grism exposure 
time of all the survey fields, and reveals many in¬ 
teresting objects. We identify more than 40 ELGs, 
including a concentration of galaxies at z ~ 0.3. 

16.2. Field J 12:19+06:49 

This field was parallel to a program aiming at a 
Virgo cluster galaxy, located in the outer regions of 
the cluster. We detect a small, low surface bright¬ 
ness spheroidal galaxy at 12:19:10.96, -1-06:47:54.0 
(J2000), which has a bright, compact source at its 
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Fig. 12.— Simulated emission line signal-to-noise ratio (S/N), as a function of galaxy magnitude (left) and 
radius (right). The solid line is the mean S/N of an [O II] emission line with an equivalent width (EW) 
of 50 A. The dashed line is same for an [O III] emission-line with EW= 150 A , corresponding to the SBl 
spectral template Kinney et al. (1996) redshifted to z = 0.7. The fraction of [O II]-sources detected at > 3cr 
(in percent) are shown in italics. 


optical center with blue-excess continuum. How¬ 
ever, it is unclear if the central object is associated 
with the galaxy. Based on the galaxy patchy faint 
background traced out to galactocentric radii of 
« 10" in both E7751T and F850LP-band images, 
the galaxy appears to be an isolated spheroidal 
galaxy, likely at the distance of the Virgo cluster. 
At the same time, we could not rule out the possi¬ 
bility that this is a Galactic source (e.g., planetary 
nebula). No record of this object exists in NED or 
SIMBAD. 

4.6.3. Field J 10:03+29:06 

This field shows a surprisingly large number 
of ELGs, despite its short total grism exposure 
time. Eleven emission line galaxies and sev¬ 
eral candidates high-redshift Lya emitters were 
observed with Keck/LRIS, as discussed in §3.2; 
however, conditions for this ground-based spec¬ 
troscopy were insufficient for obtaining useful data 
on any but the brightest galaxies. 

4 . 6 . 4 . Field J 08:06+06:43 

The ACS data for this field is part of the GO 
program 9405 (P.I. Andrew Fruchter), which tar¬ 
geted the host galaxy of supernova SN 2002LT, 


associated with gamma-ray burst GRB021211 
(Crew et al. 2003). The ACS grism spectrum 
confirms the presence of an emission line at 
~ 7450 A at the position of the probable host 
galaxy (J 08:08:59.828, -1-06:43:37.52), as previ¬ 
ously detected in a VLT/FORS2 spectrum by 
Vreeswijk et al. (2003). This line is likely [O II] 
A3727 at z=1.006. The deep multi-epoch slitless 
and direct data allow us to detect 88 ELGs in this 
field. 

5. Discussion 

We present basic data derived from the ACS 
Grism Parallel Survey. The G800L grism on ACS 
provides a unique opportunity to survey large vol¬ 
umes of the universe for faint emission lines and 
at high angular resolution. The small pixel scale 
of the ACS images and our custom software for 
deblending object spectra provides us with the 
unique opportunity to identify faint star-forming 
regions across vast cosmic epochs. The large 
dataset afforded by our “random” parallel obser¬ 
vations allows us to collect data from much larger 
area than would be possible in a single GO pro¬ 
gram. 

The ACS Grism Parallel Survey complements 
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previous and ongoing surveys of [O II], [O III] and 
Ha ELGs. Our faint flux limits allow us to probe 
deeper into the ELG luminosity function. In this 
paper we present our methodology for data anal¬ 
ysis and the first results from our survey. Our ini¬ 
tial survey of 121 square arcminutes detects 601 
emission-line galaxies at redshifts 2 < 1.6. The 
line luminosities, equivalent widths and continuum 
magnitudes suggest that we are seeing galaxies 
with a broad range of star formation rates, from 
quiescently star-forming galaxies at low redshifts 
to bright starburst galaxies at z > 1. 

The survey is biased towards compact objects 
with strong emission lines. Such galaxies tend 
to be either starbursts and/or AGNs. Follow-up 
high-resolution spectra are necessary in order to 
classify each ELG by its activity. To date we have 
obtained follow-up spectra for 11 candidates, all 
of them found to be starforming ELGs with seven 
being starbursts (see §3.2). While this subsample 
is small, we infer that the fraction of AGN in our 
sample is small and that the inferred emission-line 
luminosities can be used to estimate the co-moving 
star formation density at z < 1.6 at a high confi¬ 
dence level. 

The authors also wish to recognize and ac¬ 
knowledge the very significant cultural role and 
reverence that the summit of Mauna Kea has al¬ 
ways had within the indigenous Hawaiian commu¬ 
nity; we are most fortunate to have the oppor¬ 
tunity to conduct observations from this moun¬ 
tain. ID acknowledges the support from NASA 
HST grant GO-9468. The work of DS was carried 
out at Jet Propulsion Laboratory, Galifornia Insti¬ 
tute of Technology, under a contract with NASA. 
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Table 1 

The ACS Grism Survey Fields 


Field 

ba 

G800L 

#Exp*’ Time‘s 

#Expb 

F814W 

Time‘s 

m(a)^ 

#Expb 

F775W 

Time'^ 

m(cr)'^ 

#Expb 

F850LP 

Time‘s 


#Expb 

F606W 

Time‘s 

m(a)^ 

Prog. 

J 01 22 - 28 24 

-83 

25 

12825 




5 

2437 

29.6 

5 

2530 

28.8 




9482 

J 01 30 - 16 04 

-76 

41 

18525 




11 

5768 

30.3 

8 

3700 

29.3 




9482 

J 02 27 - 40 55 

-66 

26 

12218 

13 

6195 

30.3 










9468 

J07 26 + 6915 

28 

22 

15526 




6 

3200 

29.8 

6 

3060 

29.1 




9482 

JOS 08 + 06 43 

20 

22 

22390 

24 

11680 

30.8 







39 

15580 

30.9 

9405 

J 1003+ 29 06 

53 

27 

13481 

17 

8750 

30.5 










9468 

J 11 29 - 14 39 

44 

26 

17398 




2 

1000 

29.2 

19 

12911 

30.0 




9482 

J 12 19+ 06 49 

68 

25 

16800 




20 

10938 

30.5 

15 

9798 

29.7 




9482 

J 13 39 + 00 08 

61 

27 

13678 




8 

4460 

29.9 

8 

3760 

29.1 




9482 

J 13 58+ 62 39 

53 

42 

24319 

18 

10226 

30.7 

4 



2 

1000 

29.0 




9468 

J 15 42 - 10 46 

34 

43 

24024 




21 

11586 

30.1 

20 

10669 

29.8 




9482 


^Galactic latitude, in degrees. 

^Number of images and spectra. 

^The total exposure time, in seconds. 

^The la sky rms magnitude (Vega system), measured in the aperture of 1 x FWHM diameter (0.1 arcsec). In each filter, the number represents the average 
value for a whole area of the combined frame. 




^SYNPHOT standard Vega magnitude 

Note. —The complete version of this table is in the electronic edition of the Journal. The printed edition 
contains only a sample. 
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Table 3 

Spectroscopic Properties of Emission-Line Objects 


Id 

z 

Flux 

Hq 

EW 

S/N 

cont. 

Flux 

H^+[oni] 

EW 

S/N 

cont. 

Flux 

[Oil] 

EW 

S/N 

cont. 

R^ 

Comment.^ 

410 

1.09 









69.9± 9.9 

251.8±35.6 

7.1 

0.3 

b 


411 

0.31 

44.9± 5.8 

141.6±18.3 

7.7 

0.3 

32.7± 8.3 

61.3±15.6 

3.9 

0.5 





b 

HII? 

412 

1.15 









92.4±30.8 

92.1±30.7 

3.0 

0.8 

c 


413 

0.69 





41.7± 5.0 

338.1±40.6 

8.3 

0.1 

8.7± 2.9 

48.6±16.2 

3.0 

0.1 

b 


414 

0.71 





31.4± 1.2 

510.8±20.1 

25.4 

0.0 





a 

BCG 

415 

0.83 





60.7± 3.5 

363.3±21.1 

17.2 

0.1 

7.4± 2.5 

34.7±11.6 

3.0 

0.3 

c 

S 

416 

0.62 





64.0±21.3 

93.6±31.2 

3.0 

0.7 

91.3±30.4 

113.6±37.9 

3.0 

0.8 

c 

BCG? 

417 

0.32 

29.0± 6.1 

88.4±18.6 

4.8 

0.3 

55.5±14.1 

109.9±28.0 

3.9 

0.5 





a 


418 

0.37 





184.9±12.4 

291.5±19.5 

15.0 

0.6 





a 


419 

0.01 

111.6± 2.8 

956.7±23.6 

40.5 

0.1 









b 

BCG? 

420 

0.29 

177.4±59.1 

253.9±84.6 

3.0 

0.6 

130.1±43.4 

132.0±44.0 

3.0 

1.0 





a 


421 

1.32 









48.0±16.0 

141.6±47.2 

3.0 

0.3 

c 


422 

0.69 





115.7± 3.4 

1606.3±47.6 

33.8 

0.1 

11.4± 1.6 

132.4±18.8 

7.0 

0.1 

b 

BCG 

423 

0.85 





293.8±21.7 

359.2±26.6 

13.5 

0.8 

44.9±15.0 

29.2± 9.7 

3.0 

1.5 

b 


424 

0.98 









40.7± 5.4 

158.4±20.9 

7.6 

0.3 

b 

BCG 


Note.— All observed fluxes are in unites of 10 ergs cm ^ s and equivalent widths are in A. Fluxes and equivalenth widths are not corrected to the rest frame and absorption. 
^Reliability of line identification: ‘a’ = ’Good’, ‘b’ = ‘OK’, ‘c’ = ‘Uncertain’ 

^Symbols are as follows: ‘BCG’ = ‘Blue Compact Galaxy’, ‘S’ = ‘Spiral’, ‘Ir’ = ‘Irregular’, ‘HII’ = ‘HII region’, ‘#c’ = ‘# components’. 

Note.— The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 
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